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ABSTRACT 

We calculate the angular cross-correlation function between background QSOs and 
foreground galaxies induced by the weak lensing effect of large-scale structures. Re- 
sults are given for arbitrary Friedmann-Lemaitre cosmologies. The non-linear growth 
of density perturbations is included. Compared to the linear growth, the non-linear 
growth increases the correlation amplitude by about an order of magnitude in an 
Einstein-de Sitter universe, and by even more for lower f2 . The dependence of the 
correlation amplitude on the cosmological parameters strongly depends on the nor- 
malization of the power spectrum. The QSO-galaxy cross-correlation function is most 
sensitive to density structures on scales in the range (1 — 10) h~ x Mpc, where the nor- 
malization of the power spectrum to the observed cluster abundance appears most 
appropriate. In that case, the correlation strength changes by less than a factor of ^ 2 
when Slo varies between 0.3 and 1, quite independent of the value of For fig <> 0.3, 
the correlation strength increases with decreasing Oo, and it scales approximately 
linearly with the Hubble constant h. 



1 INTRODUCTION 

The existence of angular cross correlations between QSOs of 
moderate or high redshift with luminous foreground mate- 
rial on angular scales of J> 10' has observationally been es- 
tablished in numerous studies. Fugmann (1990) correlated 
optically bright, radio-loud QSOs with Lick galaxies and 
found a significant overdensity of galaxies around the QSOs 
of some of his QSO subsamples. Bartelmann & Schneider 
(1993b) repeated Fugmann's analysis with a well-defined 
sample of background QSOs, confirming the correlation at 
the 98% confidence level for high-redshift, optically bright 
radio-QSOs. With a similar correlation technique, correla- 
tions between the 1-Jansky QSO sample and IRAS galax- 
ies (Bartelmann & Schneider 1994) and diffuse X-ray emis- 
sion (Bartelmann, Schneider, & Hasinger 1994) were inves- 
tigated, yielding qualitatively the same results. Hutchings 
(1995) found evidence for an excess at the J> 5<r level of 
faint galaxies in fields of ~ 1' radius around seven QSOs 
with z = 2.3. 

Rodrigues- Williams & Hogan (1994) found a highly sig- 
nificant correlation between optically selected, high-redshift 
QSOs and Zwicky clusters. They discuss lensing as the most 
probable origin of the correlations, although simple mass 
models for the clusters yield lower magnifications than re- 
quired to explain the significance of the effect. Seitz & 
Schneider (1995) repeated Rodrigues- Williams & Hogan's 
analysis with the 1-Jansky sample of QSOs. They found 
agreement with the previous result for intermediate-redshift 
(z ~ 1) QSOs, but failed to detect significant correlations for 
higher-redshift sources. A variability-selected QSO sample 
was correlated with Zwicky clusters by Rodrigues- Williams 



& Hawkins (1995). They detected a significant correlation 
between distant QSOs and foreground Zwicky clusters, and 
interpreted it in terms of gravitational lensing. Again, the 
implied average QSO magnification is substantially larger 
than that inferred from simple lens models for clusters with 
typical velocity dispersions. Wu & Han (1995) searched for 
associations between distant 1-Jansky and 2-Jansky QSOs 
and foreground Abell clusters. They found no correlations 
with the 1-Jansky sources, and a marginally significant cor- 
relation with 2-Jansky sources. They argue that lensing by 
individual clusters is insufficient for typical cluster velocity 
dispersions, and that lensing by large-scale structures pro- 
vides a viable explanation. 

Bemtez & Martinez-Gonzalez (1995, 1997) found an ex- 
cess of red galaxies from the APM catalog with moderate- 
redshift 1-Jansky QSOs on angular scales < 5' at the 99.1% 
significance level. Their colour selection ensures that the 
galaxies are most likely well in the foreground of the QSOs. 
Bartsch, Schneider, & Bartelmann (1997) re-analyzed the 
correlation between IRAS galaxies and 1-Jansky QSOs us- 
ing a more advanced statistical technique which can be op- 
timized to the correlation function expected from lensing 
by large-scale structures. In agreement with Bartelmann 
& Schneider (1994), they found significant correlations be- 
tween the QSOs and the IRAS galaxies on angular scales of 
~ 5', but the correlation amplitude is higher than expected 
from weak lensing by large-scale structures. 

All these results indicate that there are correlations 
between background QSOs and foreground "light", with 
"light" either in the optical, the infrared, or the (soft) X-ray 
wavebands. Since the foreground light emission is separated 
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from the QSOs by typically a few hundred Mpc's, physi- 
cal associations are ruled out. Gravitational lensing is then 
the most likely interpretation. By the magnification bias, 
objects which are magnified by lensing are preferentially in- 
cluded in flux-limited samples. A higher than average frac- 
tion of these objects is therefore found in the background 
of matter overdensities. If light traces the overdensities, a 
correlation between foreground light and background QSOs 
can be established. 

The angular scale of the correlations of a few arc 
minutes is compatible with that expected from lensing by 
large-scale structures, and the amplitude is either consis- 
tent with that explanation or somewhat larger (Bartelmann 
& Schneider 1993a). Wu & Fang (1996) discussed whether 
the autocorrelation of clusters modeled as singular isother- 
mal spheres can produce sufficient magnification to explain 
this result. They found that this is not the case, and argue 
that large-scale structures must contribute substantially. In 
any case, the correlation scale is almost two orders of mag- 
nitude larger than the typical Einstein radius of individual 
galaxies. Hence, correlations on such scales cannot be at- 
tributed to the lensing effects of individual galaxies alone. 

If lensing is indeed responsible for the correlations de- 
tected, other signatures of lensing should be found in the 
vicinity of such sources which are correlated with foreground 
light. Indeed, Fort et al. (1996) searched for the shear ex- 
pected from weak lensing in the fields of five luminous QSOs 
and found coherent shear fields in all of them. In addition, 
they detected galaxy groups in three of their fields. Ear- 
lier, Bonnet et al. (1993) found evidence for coherent weak 
shear in the field of the potentially multiply-imaged QSO 
2345+007, which was later identified with a distant cluster 
(Mellier et al. 1994; Fischer et al. 1994; Pello et al. 1996). 

Bartelmann (1995) has analytically calculated the an- 
gular cross-correlation function between QSOs and galaxies 
assuming that (i) lensing effects are weak and (ii) the bias- 
ing hypothesis of galaxy formation (e.g. Kaiser 1984; Dekel 
& Rees 1987) holds. This study was restricted to (iii) linear 
growth of density perturbations and (iv) Einstein-de Sitter 
(EdS) model universes. The goal of this paper is to keep as- 
sumptions (i) and (ii), but give up the restrictions (iii) and 
(iv). This will enable us to study the influence of cosmolog- 
ical parameters on the correlation function, and to assess 
the importance of non-linear growth of the density pertur- 
bations. Section [| reviews the formalism, part of which is 
derived in appendix We present the results in Sect. ^| and 
discuss them in SectTTi 



2 FORMALISM 

As shown by Bartelmann (1995), the cross-correlation func- 
tion £qg(0) between high-redshift QSOs and low-redshift 
galaxies induced by weak gravitational lensing can be writ- 
ten 



£qg(4>) = ('- 1)&^W0 



(i) 



As mentioned in the introduction, the assumptions under- 
lying eq. ([]]) are that (i) lensing effects are weak (more pre- 
cisely, magnifications |/_t — 1| <C 1), and that (ii) galaxies are 
linearly biased relative to the dark matter by a factor b. s is 
the double-logarithmic slope of the differential QSO number 



counts uq as a function of flux S, nq oc S~ s , and £ M< 5 is the 
angular cross-correlation function between magnification yu 
and density contrast 8. The factor (s — 1) quantifies the mag- 
nification bias in eq. ([l]). In order to evaluate ([!]), we thus 
have to compute ^s- 



2.1 Cross-correlation of magnification and density 
contrast 

A thin light bundle is sheared by the weak gravitational 
lensing effect of the cosmic matte r distribution by the (2 x 
2) tensor -yy derived in eq. (A1C). Since lensing conserves 
surface brightness, the light bundle's magnification is given 
by the inverse determinant det -1 (Sij + jij). Dealing with 
weak lensing effects, we can assume <C 1. Therefore, 
the magnification fi = 1 + Sfi is 



1 + 5fi = [det(Sij + 7y)] 1 ?a 1 + tr(7;j) 



(2) 



Via the trace of 7y, the magnification fluctuation con- 
tains the trace of the (2 x 2) Hessian of the Newtonian 
potential fluctuations $. This can be augmented by $,33 be- 
cause the potential derivative along the light ray does not 
contribute to its deflection as long as the potential can be 
considered static while it is being passed by the light ray. 
Thus, V 2 $ is introduced, which is related to the density 
contrast 5 through Poisson's equation. In comoving coordi- 
nates, 



V 2 $ 



3»p 
2a 



(3) 



where lengths are measured in units of the Hub ble l ength 



H c, and the potential is scaled by c . With (A1C), the 
magnification fluctuation Sfi becomes 

/ fK{w')j K {w - w') 6(w') 



5/j,(w) = 3f2o / dw 



fn(w) 



(4) 



It depends on the (comoving) source distance w and on 
the direction {8-1,62) into which the light ray starts off 
at the observer. This is expressed in (kh by the vector 
= [/k(w)0i, /jr(w)02, w\j with the curvature-dependent 



radial distance /x(w) of (AS 



The magnification fluctuation as a function of distance, 
8fi(w), has to be integrated along the line of sight, weighted 
by the normalized distance distribution of the sources. Let 
Wq(w) be that distribution, then 



8p,(e 1 ,e 2 )= dw' w Q (w') Sfi(w) 

Jo 



(•») 



The upper bound on the integral is to be interpreted as w at 
z = 00, which can be finite for some combinations of cosmo- 
logical parameters. By changing the order of integrations, 
(n) can be written 



<Sa*(0i,02)= 3 fi / dwG Q (w) 



S(w) 



(6) 



with the modified weight function 

GqM = /*(«>) I™ dw' fK l W ' ~ W) W Q ( W ') . (7) 

Likewise, the density contrast must be integrated along 
the line-of-sight, weighted by the distance distribution of the 
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galaxies Wg(w) in order to get the projected galaxy num- 
ber density on the sky, according to the biasing hypothesis. 
Hence, 



S(0i 



dw5(w) Wg(w) 



(8) 



For the shape of the weight functions Wg{w) and Gq(w), 
see Fig. ^. 

Both 5 and Sp, are projections of the three-dimensional 
density contrast. The density contrast is supposed to be a 
homogeneous and isotropic random field characterized by a 
power spectrum Ps{k). Typical scales in Ps(k) are on the 
order of several ten Mpc's. The projections along the line- 
of-sight have length scales comparable to the Hubble ra- 
dius. Under these circumstances, the Fourier-space analog 
of Limber's equation can be applied (Kaiser 1992). It as- 
serts that the power spectrum P p (k) of a projection p of a 
three-dimensional field / with weight function q is related 
to the power spectrum Pf(k) of / through 



P p (k) = 



q 2 (w) 



dw -^rpr Pf 



(9) 



Note that k is a two-dimensional wave vector, while k is 
three-dimensional. The cross power spectrum P A ,a(«;) of the 
density contrast and the magnification fluctuation is there- 
fore given by 



P^,s{k) — 3fi 



dw G Q HWoH Ps 



(10) 



The cross correlation function £ M <s(^>) is the Fourier 
transform of P M a(/«), 

f°° KdK f 27T 

t»s(<f>)= / -i—^P^ii^) I d-d exp(-ire0costf) , (11) 
Jo ( 27r ) Jo 

where t? is the angle enclosed by k and <j>. The ^-integration 
yields a zeroth-order Bessel function Jo, and inserting ( |lo| ) 
leads to 



3n 



dw 



Gq(w)W g (w) 
a(w) 



^P s (k,w)Jo[fK(w)kq>] 



(12) 



It should be noted that eq. ( |l2[ ) automatically accounts 
for two effects: First, matter inhomogeneities "behind" the 
QSOs do not contribute to the gravitational lensing effects. 
This can be seen by inserting Wq = S(w — wq) into eq. (Q), 
upon which Gq oc H(wq —w), with the Heaviside step func- 
tion H(a;). Inserting this tog ether with Wg = S(w — wg) 
with wg > W Q into eq. ( p2[ ) yields = 0, as it should. 
Second, the autocorrelation of galaxies is included through 
the matter power spectrum P$(k). 

The matter power spectrum P$(k, w) depends on time, 
which is expressed in (jl2|) as a dependence on comoving dis- 
tance w. Assuming linear growth of the density fluctuations, 
P$(k,w) = Pg(k)D(w), with the power spectrum Pg(k) at 
the present epoch and a cosmology-dependent growth func- 
tion D(w). An accurate fit formula for D(w) is given in Car- 
roll, Press, & Turner (1992). In that case, the integrations 
over w and k can be separated. This yields the instructive 
equation 
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Figure 1. Examples for the galaxy weight function Wg(w) (left 
panel), and for the QSO weight functions Gq(w) adopted from 
Pei (1994; right panel). All curves are shown for an Einstein-de 
Sitter model universe, mo is the detection threshold magnitude 
for the galaxies, and zo is the lower redshift cutoff of the QSO 
sample. Line types distinguish results for different mo or zo, as 
indicated. 



kdk 
2tt 



Ps(k,w)F(k, l 



(13) 



which shows that the cross correlation function between 
magnification and density contrast is a convolution of 
the matter power spectrum with a filter function F(k,<f>), 
weighted by the wave number k. For the Einstein-de Sit- 
ter case, F(k,cj>) is plotted in Bartelmann (1995; Fig. 2). 
For small angles cf>, F(k,<f>) filters out scales smaller than 
~ 1 Mpc. Here, we evaluate eq. (|l^) numerically, including 
the non- linear growth of Pg(k,w). 

2.2 Choice of the weight functions 

We have to specify the QSO- and galaxy weight functions 
Gq and Wg- For the QSO weight function Wq, we choose a 
QSO redshift distribution of the form given by Pei (1994), 
with parameters adapted to observations in the same paper. 
With (^), this yields QSO weight functions Ga(w) for which 
examples are given in the right panel of Fig. Py. Gq depends 
on the lower redshift cutoff zo imposed on the QSO sample. 

For the galaxy redshift distribution underlying the 
galaxy weight function Wg(w), we assume a constant co- 
moving number density of galaxies. The galaxies are dis- 
tributed in luminosity according to a Schechter luminos- 
ity function with parameters v = —0.81 and L* = 1.1 x 
10 10 LoyhT 2 . We further choose a threshold apparent magni- 
tude mo for the galaxies, which we calculate from L and w 
using a if-correction derived from an effective spectral en- 
ergy distribution with index a — —1.5. Examples for the re- 
sulting weight functions are given in the left panel of Fig. [I]. 
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Our results for the cross-correlation function ^ fl s(4>) show 
that £116 (4>) depends only very weakly on the exact shape 
of the weight functions Wq(w) and Gq(w), once mo and zo 
are fixed. 



2.3 Matter power spectra 

We use in the following three different types of matter power 
spectra P$(k), viz. cold dark matter (CDM), hot dark matter 
(HDM), and mixed dark matter (MDM). For each of these 
dark-matter types, several fitting formulae have been given 
in the literature. Comparing results for the QSO-galaxy 
cross-correlation function, we found that different formulae 
for the same dark-matter type yield only marginally different 
correlation functions. We can therefore restrict the study to 
one formula for each dark-matter type only. For definiteness, 
we use those given by Efstathiou, Bond, & White (1992) for 
CDM, Bardeen et al. (1986) for HDM, and Holtzman (1989) 
for MDM. 

Semi-analytic prescriptions for the non-linear evolution 
of the power spectra have been given by Jain, Mo, & White 
(1996) and Peacock & Dodds (1996). Both rest on the con- 
jecture by Hamilton et al. (1991) that the non-linear power 
on a given scale fc n i is identified with the linear power on 
a larger scale fei which is some universal function of k n i- 
This universal function has to be identified through nu- 
merical simulations. The prescription by Peacock & Dodds 
(1996) works as long as the local double- logarithmic slope 
of the power spectrum at the scale of interest is not steeper 
than that of the CDM spectrum. This is not the case at 
small scales for an HDM power spectrum. For that, we 
rather use the evolution calculated in the framework of 
the Zel'dovich approximation by Schneider & Bartelmann 
(1995). For MDM, we use linear evolution only. To give an 
example, Fig. |i| shows the non-linear evolution of a CDM 
spectrum in an Einstein-de Sitter universe. 

One further point in specifying the power spectra is to 
normalize the shape predicted for the chosen matter model 
with measurements. There are three basic possibilities to do 
this, (i) Determine the variance of galaxy counts on a fixed 
scale, commonly 8/i _1 Mpc, and assume that light follows 
mass, with some linear bias factor b. (ii) Measure the abun- 
dance of galaxy clusters and fix the amplitude of the power 
spectrum such that this abundance is reproduced, (iii) Fix 
the amplitude of the power spectrum such that the ampli- 
tude of the cosmic microwave background anisotropies are 
reproduced. It is well known that these methods yield differ- 
ent amplitudes for the CDM power spectrum. It is possible 
to match (i) and (ii) by choosing the bias factor suitably. To 
also match (i) or (ii) with (iii) , the shape of the power spec- 
tra needs to be changed on small scales. We emphasize here 
that the angular cross-correlation function between QSOs 
and galaxies provides in principle a more direct measure- 
ment of the matter power spectrum than counting galax- 
ies or clusters, for it is based on gravitational lensing. We 
shall later demonstrate the huge effect which the non-linear 
growth of density fluctuations has on the amplitude of the 
cross-correlation function on small angular scales (Fig. fy. 

We show results for all three normalizations of the 
power spectrum, not accounting for the further constraints 
that are derived from galaxy counts, galaxy autocorrelation 




0.001 0.01 0.1 



k [h/Mpc] 

Figure 2. Comparison between the linear (dotted curves) and 
non-linear (solid curves) evolution of a CDM power spectrum, 
according to the prescription given by Peacock &: Dodds (1996) 
for an Einstein-de Sitter model universe. Non-linear power builds 
up on small scales at the expense of power on intermediate scales. 
The linear spectra are normalized to <rg = 1. 

functions, or the cluster abundance. In that way, our results 
can easily be scaled to the desired normalization. 



3 RESULTS 

As shown by Bartelmann (1995) for the special case of a 
linearly evolving density field in an Einstein-de Sitter uni- 
verse, the shape of the cross-correlation function ^s(4>) de- 
pends sensitively on the type of the dark-matter particles. 
This property of £ M a(0) continues to hold for all the gener- 
alizations we calculate here. We use the results obtained for 
linear density evolution in an Einstein-de Sitter universe as 
a reference for the influence of cosmological parameters and 
non-linear density evolution. 

3.1 Dependence on type and evolution of the 
dark matter 

Figure [| demonstrates the different shapes of ;((p) for dif- 
ferent matter power spectra. In agreement with the ana- 
lytic description of Bartelmann (1995), the cross-correlation 
function (left panel of Fig. ^) peaks for — > 0, whereas 
it flattens off at small 4> for HDM. Power-series expansion 
yields oc <j> + 0(</> 2 ) for CDM, and ^ s <x cp 2 + O^ 4 ) for 
HDM (Bartelmann 1995). The reason for this is the different 
asymptotic behavior of Ps(k) for k — > 00: the cut-off in the 
HDM power spectrum is responsible for the flattening of the 
HDM cross-correlation function. As expected, £ M a for MDM 
falls between these two cases, in that it flattens at smaller 
than for HDM. 
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Figure 3. The cross-correlation function £uc5(^0 is shown here 
for an Einstein-de Sitter universe. Left panel: results for different 
types of linearly evolving dark-matter power spectra (CDM, solid 
curve; MDM, dotted curve; HDM, dashed curve). All spectra were 
normalized to as = 1. While the cross-correlation function peaks 
towards <f> — ► for CDM, it flattens for HDM and MDM due to 
the cut-off at large wavenumbers k in those spectra. At large an- 
gles (j>, the HDM and MDM correlation functions lie above the 
CDM results. Right panel: the effect of the non-linear evolution 
of the CDM power spectrum on (</>), for three different nor- 
malizations, as indicated by line type. The COBE normalization 
yields the largest amplitude, the cluster normalization the small- 
est. The linear results are shown as the thin lines for comparison. 
For this choice of the cosmological model, non-linear evolution 
increases the correlation amplitude £ M a(0) by about an order of 
magnitude, quite independent of the normalization, mp = 20.5 
and zo = 0.3 were chosen. 



The cross-correlation function for the different normal- 
izations of the CDM power spectrum is shown in the right 
panel of Fig. |j. Clearly, the high normalization to the COBE 
results yields the strongest correlation, while the normaliza- 
tion to cluster counts yields the weakest. For small angles cf>, 
the window function F(k, <f>) (see Bartelmann 1995 for de- 
tails) drops to zero at the scale of galaxy clusters and smaller 
structures. Hence, the integral (h2|) yields the strongest cor- 
relation for the COBE normalization, because the COBE 
normalization leads to much more power on cluster scales 
than the cluster normalization. 

The non-linear growth of the density fluctuations also 
most strongly affects the power on scales of clusters and 
smaller (see Fig. ^). Because of that, the influence of non- 
linear density evolution on the cross-correlation function is 
huge. The right panel of Fig. ^contrasts the cross-correlation 
function for non-linear and linear density evolution of the 
power spectra. For the CDM model in an Einstein-de Sit- 
ter universe, the amplitude, f M j(0), of the cross-correlation 
function increases by about one order in magnitude, quite 
independent of the normalization. Depending on the nor- 



malization, the increase can be even larger in other model 
universes. Since ^^((p) remains unaffected on large angular 
scales, the peak for (j> — > is sharpened by non-linear evo- 
lution, in tentative agreement with observations (Bemtez & 
Martinez-Gonzalez 1995; 1997). Non- linear evolution of the 
power spectra is, therefore, crucial to match theoretical ex- 
pectations and measurements. 

Non-linear evolution also increases the HDM correlation 
amplitude. However, the flattening of the HDM correlation 
function for <j!> — > 0, and its much lower amplitude than for 
CDM, appears to be ruled out by observations, and we will 
therefore concentrate on CDM spectra in the following. 

3.2 Dependence on cosmological parameters 

The change of the cross-correlation function with different 
cosmological parameters shows a huge dependence on the 
normalization of the power spectrum. This is due to the fact 
that the three normalizations themselves depend differently 
on the cosmological parameters. 

To understand the dependence of the correlation func- 
tion on the cosmological parameters a bit better, we adapt 
an analytical expression by Bartelmann (1995) to our dis- 
cussion. Using a CDM model spectrum, and assuming linear 
density evolution, Bartelmann (1995) obtains 

^(0) = ^lMxo (14) 

for an Einstein-de Sitter universe. Here, Acdm is the ampli- 
tude of the power spectrum, ko is the wave number of the 
peak in the spectrum, and Kq is an integral over the weight 
functions introduced in Sect. 2.2. In non-EdS cosmologies, 
we get an additional factor Qo from Poisson's equation, and 
fco oc Qo h in units of the inverse Hubble length. Using this, 
the dependence of the amplitude £ M a(0) (for linear evolution) 
can be written as 

£ M a(0) oc Q(5 h? x [Acdm] x [K ] , (15) 

where the square brackets around a symbol denote the de- 
pendence of the corresponding quantity on the cosmological 
parameters. Fig. ^ illustrates the dependence of £^(0) on h 
and ilo for various cosmological models and normalizations 
of the power spectrum. 

We looked in detail at three different cosmological mod- 
els. These are: an Einstein-de Sitter universe with two differ- 
ent values of the Hubble constant Ho, an open low-density 
model with Slo = 0.3, and S7a = 0, and a spatially flat low- 
density model with Slo = 0.3 and Q,a ~ 0.7. The correlation 
functions £ M a (<f>) for these cosmologies and different normal- 
izations are shown in Fig. |B|. The left panel of Fig. ^ illus- 
trates the dependence on Qq of the correlation amplitude 
Cqg(0) for cluster normalization. We restrict the discussion 
of the dependence on cosmological parameters to non-linear 
models only. 

In an Einstein-de Sitter universe, the correlation is 
strongest for the COBE normalization. This changes in low- 
density universes, where we find a stronger correlation for 
the cluster- and the as = 1 normalizations. Generally, the 
correlation amplitude increases with Qo , except for the clus- 
ter normalization. There, the correlation amplitude changes 
only mildly with Qo as long as Qo iS 0.3, and it increases 
when Qo decreases even further. Likewise, a larger Hubble 
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Figure 4. The dependence of the correlation amplitude ^ M ^(0) 
on the Hubble constant h and the matter density Qq is shown for 
different cosmological models and different normalizations of the 
CDM power spectrum. Non-linear evolution was assumed. Left 
panels: dependence on h; right panels: dependence on Ho; upper 
panels: open models (Q\ = 0); lower panels: spatially flat models 
(Oo + = 1)- Different normalizations are distinguished by line 
type, as indicated. The COBE normalization yields the strongest 
dependence on both h and Qo- Generally, § M ,5(0) decreases with 
decreasing h and decreasing Qo, but note that the cluster normal- 
ization reverses that trend, especially for the open models (upper 
right panel). All curves were calculated assuming mo = 20.5 and 
z = 0.3. 

constant yields a higher correlation amplitude. Both these 
effects are mainly due to the fact that the peak in the power 
spectrum shifts towards larger wave numbers with increas- 
ing flo and Ho, and the additional factor S7q from Poisson's 
equation renders a stronger dependence on flo than on Hq. 
Depending itself quite strongly on f2o, the cluster normal- 
ization can reverse the general trend with Qq. 

Approximate relations between the amplitude £^(0) 
and S7o and Ho are shown in Tab. |l| They are in goo d quali- 
tative agreement with the expectations from eq. ( |15[ ) . Figure 
^ and Tab. [j] also show that the strongest dependence on the 
cosmological parameters is obtained for the COBE normal- 
ization. The dependence on the cosmological constant S2a 
for fixed Qo is much weaker than for the other parameters, 
which reflects the fact that $2a affects only the space-time 
curvature, which enters through Ko in eq. ( |l5| ) . 

3.3 Other model parameters 

Quite generally, the dependence of the correlation function 
on all the other parameters entering our calculation is weak. 
As mentioned in Sect. 2.3, there is no significant change of 
the correlation function upon using different spectra for the 
same dark-matter model. Small changes in the amplitude 
of the angular cross-correlation function result from varying 



Figure 5. The angular cross-correlation function s (e>) is plot- 
ted for different values of the cosmological parameters. In all pan- 
els, curves for the three normalizations are shown, distinguished 
by line type as indicated. Top left panel: spatially fiat low-density 
universe with h = 0.5, Qo = 0.3, and Qa = 0.7; top right panel: 
open low-density universe with h = 0.5, Qo = 0.3, and ^a = 0; 
bottom panels: Einstein-de Sitter universes with h = 0.5 (left) 
and h = 0.7 (right). All curves are for CDM and for mo = 20.5 
and zo = 0.3. 

the parameters zo and mo of the weight functions introduced 
in Sect. 2.2. This is due to the fact that varying the mean 
redshift of the galaxies (through mo), or the QSO redshift 
distribution (through zo) changes the lensing efficiency. Ex- 
amples are shown in the right panels of Fig. ^j. The upper 
right panel illustrates how the correlation amplitude changes 
if the mean redshift of the galaxies is shifted beyond the 
maximum in lensing efficiency for a fixed QSO redshift dis- 
tribution. The lower right panel shows how the correlation 
amplitude increases with increasing QSO redshift. 

The fact that ^ /J i(<^) peaks for moderately faint galaxies 
rather than for the faintest can easily be understood. The 
galaxies are simply tracers of the density inhomogeneities 
that act as lenses on the QSO samples. At fixed QSO red- 
shift, there is an intermediate redshift at which the lenses 
are most efficient. This redshift is somewhat reduced by the 
redshift dependence of structure growth. Allowing for fainter 
galaxies selects more distant density inhomogeneities which 
are less efficient lenses, making the cross-correlation function 
^5 drop. 



4 SUMMARY AND DISCUSSION 

We have extended an earlier study of the angular cross- 
correlation function £qg(</>) between QSOs and galaxies 
(Bartelmann 1995) in three respects. First, we generalize 
from the Einstein-de Sitter to arbitrary Friedmann-Lemaitre 
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Table 1. This table summarizes the dependence of the corre- 
lation amplitude (0) on the cosmological parameters h and 
Qrj; for non-linearly evolved CDM spectra. We fitted parabolae 
in ln£ M ,5 (0)-ln Q space to the curves in Fig. W, with Q either h or 
Q - Hence, we fit ln^(0) = A (In Q) 2 + B lnQ + C, and give the 
coefficients A, B and C in the table. When h is varied, we take 
Qo = 0.3, Ha = for the open models, and Qq = 0.3, Ha = 0.7 
for the spatially flat models. When f2o is varied, we take h = 0.5. 
The dependences differ with the normalization, hence we give re- 
sults for the three normalizations we have chosen. For the cluster 
normalization, we also give results for the correlation amplitude 
Cqg(0), which differs from §^(0) by the bias factor b = a^ 1 . 
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cosmological models, second, we allow for the non-linear evo- 
lution of the density fluctuations, and third, we construct 
more realistic redshift distributions for galaxies and QSOs. 
The QSO-galaxy correlation function is expressed in terms 
of the correlation function between magnification and den- 
sity contrast, £ M 5 (</>)• 

We confirm that the shape of the cross-correlation func- 
tion is very sensitive to the type of the dark-matter power 
spectrum: While it peaks for (f> — > in the case of CDM, it 
flattens for HDM and MDM. Measurements of £qgO) indi- 
cate that the correlation peaks for </!>—> 0, in contradiction 
to the flattening expected for HDM. Given further the much 
lower correlation amplitude for HDM compared to CDM, we 
focus on the CDM model. 

We summarize our main results as follows: 

(i) Non-linear evolution of the density fluctuations leads 
to an increase in the correlation amplitude £^(0) relative to 
the linear case by about an order of magnitude, or more. 

(ii) ^iisi'P) depends on the cosmological parameters Qo, 
Ho, and Qa in different ways: 

• Qo enters through Poisson's equation, through the 
normalization, the shape, and the evolution of the power 
spectrum, and through space-time curvature. The loca- 
tion of the peak in the power spectrum is proportional 
to Qo h = r. The larger F is, the larger is the overlap 
between the power spectrum and the window function 
that filters out the scales relevant for £ M a. Therefore, the 
correlation generally becomes stronger with increasing 
Qo, except for the cluster normalization. 




0.04 - 



Figure 6. Left panel: amplitude of the angular cross-correlation 
function between QSOs and galaxies, £qg(0), as a function of Qo. 
The curves are for cluster-normalized, non-linearly evolving CDM 
in open and spatially flat model universes. Note the difference 
to the correlation amplitude £ M a(0) between magnification and 
density contrast: We have assumed a bias factor of b = er^ 1 here. 
Upper right panel: dependence of §Qg(0) on the threshold galaxy 
magnitude mo; lower right panel: dependence of £qg(0) on the 
minimum QSO redshift zq- While §qg(0) increases monotonically 
with zo, it reaches a maximum for mo ~ 19 and decreases for 
fainter galaxies. The curves in the right panels are calculated in 
an Einstein-de Sitter universe with cluster normalization of the 
power spectrum, zo = 0.3 in the upper right panel and mo = 20.5 
in the lower right panel. 



• The Hubble constant enters through the shape pa- 
rameter of the power spectrum, F = Qo h. As for Qo, the 
correlation strengthens with increasing Ho- 
rn The cosmological constant Qa enters only weakly 
through space-time curvature and the evolution of the 
power spectrum. 

(iii) The method by which the power spectrum is nor- 
malized plays a critical role, both for the amplitude of 
the correlation and its dependence on cosmological param- 
eters. When the spectrum is normalized to the COBE mea- 
surements, the dependence on cosmological parameters is 
strongest. For the cluster normalization, the dependence on 
Qo is weakest, and it can even be reversed for small Qo, 
because the general decrease is compensated by the Qo de- 
pendence of the normalization. 

(iv) Other parameters entering the calculation of the cor- 
relation function have a comparably modest influence. Fac- 
tors of 1 .5 — 2 can be gained or lost by varying the galaxy de- 
tection threshold mo or the minimum QSO redshift zo- The 
exact shape of the galaxy- and QSO redshift distributions 
is largely irrelevant, and different representations published 
in the literature for the same type of power spectrum yield 
the same correlation functions. 
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fail to be applicable when /i becomes noticeably larger than 
unity, /i J; 1.5, say. This is the case for QSOs closer than ~ 3 
Einstein radii to a cluster core. Depending on cosmological 
parameters, QSO and cluster redshifts, and on the velocity 
dispersion of the cluster, ~ 3 Einstein radii correspond to 
~ 1 — 2'. That is to say that the expected QSO-galaxy cor- 
relation function calculated here should fall below the true 
correlation at <j> 1 — 2'. 

The most accurate measurement of £,qg(4>) has so 
far been obtained by Bem'tez & Martinez-Gonzalez (1995, 
1997). They find a significant correlation at angular scales 
(j> < 5', and a result compatible with zero at larger <j>. We 
can therefore not yet compare our results to the innermost 
points of the observational data. This must be postponed 
until £qg(<^) has significantly been detected at angles <; 3'. 
We note that especially the Sloan Digital Sky Survey will 
provide an ideal data base for this kind of analysis. 

A comparable study was performed by Sanz, Martinez- 
Gonzalez, & Bem'tez (1997) while the work on this paper was 
completed. In a forthcoming paper, these authors will per- 
form a detailed comparison of the data and the theoretical 
results. 



Figure 7. The range of correlation functions ^qg(</') that can 
be covered by varying cosmological parameters while keeping the 
normalization of the power spectrum fixed to the cluster abun- 
dance. As Tab. [l] and Fig. Q show for that normalization, (<p) 
changes by less than a factor 2 for Qo ft 0.3, and it scales with 
Hq approximately linearly. Note that we are plotting §QG rather 
than here. We have assumed s = 2.8 (cf. Bartelmann 1995) 
and b = cr^ 1 to relate the two; see eq. (Q). The cluster normal- 
ization requires roughly as ~ 0.6 Og' 6 , hence b oc Hq -0 ' 6 . See also 
the curves in the left panel of Fig. ^j. We have further assumed 
mo = 20.5 and zo = 0.3, corresponding to the values chosen in 
Bem'tez & Martinez-Gonzalez (1995, 1997), whose data points are 
added for illustration. 
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We have argued that, for small angular scales (j>, ^ns(<j>) 
is most sensitive to the power on cluster scales. This is be- 
cause the filter function F(k, <f>) introduced in eq. (^) drops 
to zero on smaller scales, and the further factor k in the inte- 
grand of eq. ( Jl3| ) suppresses the largest scales. Lacking more 
precise information on the true shape and normalization of 
the matter power spectrum on the scales relevant for the 
QSO-galaxy cross-correlation function, the normalization to 
the observed cluster abundance appears most appropriate 
for our purposes. The dependence of £ M a(0) especially on the 
density parameter is then fairly weak as long as flo <; 0.3. 
As Tab. jl] and Fig. ^jshow, £^(0) changes by at most a fac- 
tor of Ss 2 across that range of ilo, and the dependence on 
Hq is approximately linear. Figure R illustrates the range of 
QSO-galaxy cross-correlation functions £qg(0) that can be 
covered by varying cosmological parameters while keeping 
the normalization of the power spectrum fixed to the cluster 
abundance. For illustration, we overplotted Fig. |^ with the 
observational data points by Bem'tez & Martinez- Gonzalez 
(1995, 1997). 

While non-linear structures on clusters scales are accu- 
rately represented by the non-linear power spectrum, our 
approach does not accurately represent the strong lensing 
that sets in close to cluster cores. Since we have linearized 
the magnification factor /i in eqs. (Q) and (^|), our results 



APPENDIX A: THE COSMOLOGICAL SHEAR 
TENSOR 

We summarize in this appendix the derivation of the shear 
tensor by which the lensing effects of cosmic material deform 
thin light bundles. We start from the conformal Friedmann- 
Lemaitre (FL) metric, 

ds 2 = a 2 iff) \dr/ 2 — dw 2 — /k(w) do; 2 ] , (Al) 

with the scale factor a, the conformal time rj, the comov- 
ing distance w, and the comoving angular-diameter distance 
/ic (w), which depends on the spatial curvature K. r\ is re- 
lated to the cosmic time t through dt = adr/ with the bound- 
ary condition rj = at t — 0. Measuring lengths in units of 
the Hubble length Hq 1 c, and setting the scale factor at the 
present epoch to unity, K can be written in the usual form 



K = Q + Q A - 1 



(A2) 



where flo is the matter-density parameter, and Qa is the 
density parameter attributed to the cosmological constant, 
both taken at the present epoch. 

Consider n ow a fiducial light ray propagating through 
the space-time (Al), starting off at the observer into direc- 
tion (#i, 82)- The affine parameter A along the fiducial ray is 
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chosen such that dA = — adt. The transverse separation vec- 
tor £ of a neighboring light ray from the fiducial ray changes 
with A according to 



dA 2 



(A3) 



(Gunn 1967; Blandford et al. 1991j_Seitz, Schneider, & 



Ehlers 1994). We now substitute in (A3) the comoving sep- 
aration vector x = a _1 £ fo r £, and the comoving distance 
w for A. From ds = and (Al), dw — — drj for radial light 
rays. Since dA = —adt = —a dr/, dw = a~ 2 dA. After some 
manipulation, these substitutions yield 



dw 2 



+ Kx = 



(A4) 



The solutions of this oscillator equation are the usual linear 
combinations of trigonometric (for K > 0), hyperbolic (for 
K < 0), or line ar f unctions (for K = 0). 

Equation (A4) describes the change in comoving trans- 
verse separation between two light rays on cosmic scales, 
~ Hq 1 c. Locally, space time is perturbed by density fluc- 
tuations with Newtonian potential $. Relative to an unper- 
turbed light ray, such perturbations deflect light according 
to 



-2Vi$(ra) 



(A5) 



d 2 x{w) 
dw 2 

(e.g. Schneider, Ehlers, & Falco 1992; Narayan & Bartel- 
mann 1997) where we have scaled $ by c 2 and used that 
locally, dw = —cdt. w is the vector [/k(w)#i, /k(w)#2, w]. 
The gradient of $ has to be taken perpendicular to the ac- 
tual light ray, but for the small deflections we consider here 
it can be taken perpendicular to the unperturbed ray. 



Equations (A4) and (A5) describe light deflection on 
completely different spatial scales. To combine the effect of 
the global curvature of space time with the light deflection 
by isolated density fluctuations, we can thus sim ply a dd the 
right-hand side of (A5) to the right-hand side of (A4). Writ- 
ten in components, the result is 



d 2 Xj(w) 
dw 2 



+ Kx,(w) = -2$,, (w) , 



(A6) 



where the index i preceded by a comma denotes the partial 
derivative with respect to Xi. Equation (A.6) is solved by 



Xi(w)=f K (w)e-2 dw'f K (v)-v!')9,i(vf) , (A7) 
Jo 

with the boundary conditions that the unperturbed fiducial 
ray and the perturbed ray both start at the observer into 
directions separated by 8, hence 5?i(0) = and dxi/dw — 8 
at w — 0. fic(w) is the comoving angular-diameter distance 
of the metric (Al), 



K~ 

w 

\K 



1/2 



sin(if 1/2 w) 



-1/2 



(K>0) 
(K = 0) 
h(\K\ 1/2 w) (K<0) 



(A8) 



Recall that (A7) is the comoving separation of a perturbed 
light ray from an unperturbed one. The net angle a by which 
the perturbed light ray is deflected is therefore simply 

it t \a t I \a i x l {w)-f K {w)8 
a i [f K {w)8 1 , fK(w)8 2 ,w] = f K ( w ) 



= -2 fd /f"^ ,,^). (A9) 

(Seljak 1996) . Consequently, the shear exerted on a thin light 
bundle is 

da t [w, f K (w)8i, f K (w)9 2 ] 



7«(«0 = 



BO, 



= 2 f\ w/ M^)M W - w ') , 

Jo fcM 

(A10) 

This result agrees with that of other recent studies which 
focus on different lensing statistics (Seljak 1996; Jain & Sel- 
jak 1996; Bernardeau, van Waerbeke, & Mellier 1996; Kaiser 
1996). 
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